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1. Introduction 

In recent years several analytical models have been developed to predict 
microwave scattering by trees and forest canopies [1*6]. These models contribute to the 
understanding of radar backscattcr over forested regions to the extent that they capture 
the basic interactions between microwave radiation and tree canopies, understories, and 
ground layers as functions of incidence angle, wavelength, and polarization. The Santa 
Barbara microwave model backscattcr model for woodland (i.e. with discontinuous tree 
canopies) combines a single-tree backscattcr model and a gap probability model [5-6], 
Comparison of model predictions with synthetic aperture radar (SAR) data at L-band (A. 
= 0.235 m) is promising [7], but much work is still needed to test the validity of model 
predictions at other wavelengths. Here we test the validity of the model predictions at 
P-band (A. = 0.68 m) for woodland stands at our Mt. Shasta test site. 

2. Study area 

Extensive ground data have been collected to support our SIR-B (Shuttle Imaging 
Radar) and SIR-C/X-SAR studies. The site is located to the southeast of Mt. Shasta, 
California (41°18’ N, 122°05’ W), and spans elevations from 1160 to 1220 m. The 
forest stands under investigation arc natural, are on level ground, and are dominated by 
pondcrosa pine (Pinus ponderosa) or pine mixed with white fir (Abies concolor). Under 
story vegetation is sparse and consists primarily of perennial grasses and forbes. The 
litter layer can reach a depth of about 0.1 m. Three ponderosa pine stands ("Sp2", "St2”, 
and "Sll 1") in the site were chosen for this study because adequate ground data for these 
stands are available to provide inputs to the model. The three stands have different jstand 
densities and trunk diameters at breast height (dbh). Sp2 has trees of 25.0 m h and 
0.77 m dbh, with stand density of 23 trees per hectare. "Sill" has trees of 27.0 m h and 
0.46 m dbh, with stand density of 228 trees per hectare. "Sl2" has trees of 25 m h and 
0.42 m dbh, with stand density of 303 trees per hectare. The soil of the stands is derived 
from recent alluvial deposits of volcanic ash, and the ground surface is smooth at P- 
band. Other ground data used as model inputs can be found in [8]. 

3. SAR data 

SAR data were acquired by the Jet Propulsion Laboratory (JPL) airborne SAR 
overflights on 6 September 1989 after a dry summer. The SAR data were processed and 
calibrated by JPL. The estimated calibration uncertainty of the backscatter is +1.0 dB 
[3], The estimated calibration uncertainly of VV-HH phase difference is ±10°, and the 
estimated calibration uncertainty of HH and VV correlation coefficient is ± 0.1 (personal 
communication, Freeman 1992). We received the standard 4-look compressed data with 
pixel spacing of 12.1 m (azimuth) and 6.7 m (slant range). For the three stands, the 
means of the backscatter (HH, HV, and VV polarizations), the VV-HH phase difference, 
and the HH and VV correlation coefficient were extracted from the SAR data. 
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4. Comparison of model results with SAR data at P-band 

Model simulations for stands "St2", "Stl 1", and "Sp2" were carried out at the same 
incidence angle (9 0 ) as the SAR data. We have 8 SAR data takes covering stand "St2" 
(0 O ranges from 22° to 54°), 5 data takes covering stand "Stl 1" (0 O from 25° to 45°), and 
5 data takes covering stand "Sp2" (0<, from 27° to 47°). The HH, HV, and VV back- 
scatter, and the HH and VV correlation coefficient are shown in Figure 1. The model 
makes good predictions of the HH backscatter for stands "St2" and "Stll” (Figure la), 
the VV backscatter for all three stands (Figure lb), and the HV backscatter for stand 
"Stl 1" (Figure 1c). The model overestimates the HH backscatter for stand "Sp2" (Figure 
la). The model underestimates the HV backscatter for stand "Sp2" and stand "St2" (Fig- 
ure lc). We attribute the underestimate to the surface model, in which we model only 
the co-polarized surface backscatter. Since the soil is dry and the radar wavelength is 
long, deep penetration into the soil is expected. Because "Sp2" has the low stand density 
and "St2" has small trees, more soil surface is exposed in these stands. Thus, the back- 
scatter from the sub-surface may be a major contributor. 

The model predicts well for the HH and VV correlation coefficients of stands 
"St2" and "Stll", but overestimates the coefficient of stand "Sp2" (Figure Id). The 
overestimate can be attributed to inadequate modeling of surface backscatter. The sur- 
face model used gives completely correlated backscatter for HH and VV because identi- 
cal surface parameters are used for all the simulated pixels in a stand, and because the 
surface model predicts zero cross-polarized backscatter. For the sparse stand ”Sp2", the 
surface backscatter is the major scattering source. Thus, strong and completely corre- 
lated surface backscatter is predicted by the model. 

The model predicts a phase difference of ~ 210.0° for the dense stands "St2" and 
"Stll". Since the model phase difference and the SAR observed phase differences are 
close to 180° for stands "St2" and "Stll", there is some indication that the VV-HH phase 
difference from double-bounce trunk-ground interactions (180°) contributes to the total 
phase difference. 

To investigate the model prediction quantitatively, we test the null hypothesis of 
U 0 : \imodei = Vsary where \imodei is the modeled mean and } x sar is the mean of the SAR 
data. At the 5% significance level, the U 0 tests of the HH, HV, and VV backscatter, the 
correlation coefficient of HH and VV polarizations, and the VV-HH phase difference are 
accepted. 

5. Conclusions 

For P-band, the model predicts well compared with SAR data for HH, HV, and VV 
backscatter, VV-HH phase difference, and HH and VV correlation coefficient for stands 
with medium-to-high density (228 - 303 trees per hectare). When stand density is low 
(23 trees per hectare), model performance becomes unacceptable: the model overesti- 
mates HH and VV backscatter, underestimates HV backscatter, and overestimates HH 
and VV correlation coefficient. 
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Figure 1 , Model results vs SAR data at P-band 



(a). Model HH backscatter 


(b). Model W backscatter 



(c). Model HV backscatter 


(d). Model correl. coef. of HH & W 


Solid lines indicate the equality of modeled means to those of SAR data. Dotted 
lines show ± 1.0 dB calibration uncertainty for the HH, HV, and VV backscatter (Figure 
la-c) and ±0.1 for the HH and VV correlation coefficient (Figure Id). and "x" 

represent data points of stands "Sp2”, "St2", and "Stl 1", respectively. 
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